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ABSTRACT 

Procedure for calculating various interactions of neutral and charged 
particle efflux from ion thrusters with spacecraft surfaces is outlined. 
Calculation details are referenced. Details of charge -exchange calcula- 
tion in the near -field of the thruster are presented. 
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PROCEDURE FOR ESTIMATING EFFECTS OF ION BEAM 


INTERACTION WITH SPACECRAFT 
by Thaine W. Reynolds 
Lewis Research Center 

SUMMARY 

Each spacecraft layout involving electrostatic thrusters presents 
its own unique set of conditions which should be systematically checked 
for possible interaction between the ion beam and spacecraft compo- 
nents. For this reason, this report presents an outline of the proce- 
dure that has been followed by the author at the Lewis Research Center 
to check possible spacecraft contamination from ion thruster operation. 

Details of calculation procedure are not presented if they are 
available in already published references. The references are cited 
in each instance. 

Details and graphical results of charge exchange calculations are 
presented since this material has not been previously published. 

INTRODUCTION 

The possible problems associated with the effects of the thruster 
exhaust on various spacecraft components for spacecraft employing 
electrostatic thrusters have been enumerated and discussed in previous 
publications (refs. 1 to 8). Since each spacecraft design or layout is 
unique, it is necessary to follow through some general procedure to 
check all interaction possibilities to be assured that no detrimental ef- 

_ _ . -fects will be encountered. - - - - - - 

The purpose of this report is to summarize a general procedure 
that has been followed at Lewis Research Center. It includes all the 
interaction effects associated with particle efflux from the thrusters 
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that have been considered as possibly detrimental. Calculation details 
are not presented if they are available in already published references. 
The appropriate references are cited. 

A procedure used in calculation of charge -exchange neutral and ion 
formation rates is discussed in some detail, since the method used has 
not been published previously. 

ANALYTICAL PROCEDURE 

The following sketch illustrates the principal types and sources of 
particle efflux from ion thrusters that are of sufficient magnitude to af- 
fect possible spacecraft components and are considered in the analysis. 



Primary 

surface 


Types of particle flux illustrated are: 

(a,) Main beam propellant ions (fast Hg + ) 

(b) Propellant neutrals resulting from beam inefficiency (slow Hg°) 
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(c) Propellant ions from charge exchange between (a) and (b) 

(slow Hg + ) 

(d) Grid metal neutrals sputtered by (c) ions (slow M°) 

(e) Grid metal ions from charge exchange between (a) -ions and 

(d) -neutrals (slow M + ) 

(f) Neutral atoms sputtered from a primary surface by (a) -ions 

(slow M°) 

The problems associated with these fluxes may arise from sputter 
ing caused by the impingement of the higher energy particles, or from 
build-up of a layer of condensate of the arriving material. 

Items (a), (b), and (d) are of concern only_with surfaces in line -of- 
sight with the thruster exit plane (referred to as primary surfaces). 
Items (c), (e), and (f) may be concerned as well with surfaces which 
are not in line -of -sight of the thruster exit plane (referred to as 
secondary surfaces). 


Information Required 

The following items of information are required: 

(1) Design layout of the spacecraft, including ion engine and all 
components that might be affected by sputtering or condensation of 
material. 

(2) Thruster data 

(a) Size (radius R Q ) 

(b) ' Operating data: total propellant flow, propellant utiliza- 

tion efficiency, accelerator impingement current, 
accel-decel voltages 

(c) Experimental ion-beam profile at operating- conditions, 

taken sufficiently far downstream to be representative 
of the far -field profile (i.e. , at Z/R Q > 8) 
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PROCEDURE 

The following list constitutes essentially the procedure followed . 
The references cited in each item contain the detailed discussions of 
that particular phase of the calculations. 

(1) The far -field experimental ion beam profile data is fitted by 
the equation 

Ke). = e -[Mi-cos e)] n ^ 

3(0) 

to determine the constants n and A (ref. 9). This equation, with the 
determined values of n and A, is then used in all subsequent calcu- 
lations requiring use of an ion beam density value at any location be- 
yond an axial distance, Z/R q > 8. Calculations requiring values of 
ion beam density at locations closer than Z/R q = 8 are treated differ- 
ently as will be noted later . . 

(2) Using equation (1), constant current density lines are plotted 
on the spacecraft layout (or on a transparent overlay to scale) to check 
ion -impingement current values on any primary surface. (Plots are 
illustrated by fig. 9 of ref. 9.) 

(3) The ion impingement currents, determined in item (2), and a 
knowledge of the sputtering coefficient (refs. 10 and 11) are combined 
to determine the sputtering rate for any primary surface. 

(4) View factors (refs. 6, 12, and 13) from primary surfaces to 
any sensitive secondary surface are determined. The product of these 
view factors and the sputtering rate determined in item (3) are checked 
to see if the arrival rate of sputtered material from the primary: sur - _ 
face is of sufficient magnitude to cause any problem (ref. 7). 

(5) View factors (refs, as item 4) from the thruster to any primary 
surface are determined. (A plot of constant view factor lines similar 
to that made for ion current density mentioned in item 2 is a convenient 
way of exhibiting these.) These view factors are applied to neutral 
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emission rates of both neutral propellant atoms and sputtered grid 
neutrals. The corresponding arrival rates are compared with desorp- 
tion rates at the surface temperatures to see if condensing conditions 
exist (ref. 6). If so, determine whether the condensation rates are 
high enough to cause a problem over the mission duration (refs. 2, 

7, and 8). 

(6) If the neutral arrival rates calculated in item (5) are not high 
enough to condense on primary surfaces, their subsequent arrival 
rates on any secondary surface are checked (ref. 7). 

(7) Charge -exchange formation rates are calculated for both pro- 
pellant atoms and sputtered grid metal atoms. The procedure and 
some graphical results are presented in the CHARGE EXCHANGE 
section of this report. 

The subsequent trajectories of these charge exchange ions is a 
matter of conjecture at this time. Since these ions are principally 
low velocity they are expected to follow whatever electric field lines 
exist in their vicinity. Several types of assumption have been made 
as to their trajectories, including: (a) that the charge exchange ions 
travel radially outward (i. e. , radially from the ion beam centerline) 
from the point of their formation, (b) that they travel normal to some 
defined "beam -spreading” angle, and (c) that they disperse spher- 
ically symmetrically from the point of their formation, essentially as 
if no general electric field existed to guide them elsewhere. Exper- 
imental data in this area is needed to better define this problem. 

CHARGE EXCHANGE 

- - This section gives the-details of the charge exchange calculation * 
for ion thrusters. 

A major portion of the total charge exchange reactions takes place 
within a few thruster radii of the thruster exit. The use of the sim- 
plified far -field relations for current density and neutral density are 
thus not justified in this calculation. In the following development ap- 
plicable near -field relations are used. 
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The following sketch shows the definitions of variables used. 


^ Profile plane 
\ 



The total number of ions exchanging charge in an elemental volume 
at (r, z) is 

dN = ct 10 ju^ dV = a 1Q j dS dZ (2) 

(Symbols are defined in the appendix. ) In the cylindrically symmetric 
case, using ring elements, dS = 2jrr dr. 

In this development, a specific distribution for the ion emission 
rate across the emitter face is assumed, namely, 


j(o,o) 



(3) 


The ion arrival rate at the point (r,z) is accordingly given by the follow- 
ing expression (developed in appendix B. of ref.. 9) . 


- 2. 16 C(n, A) n , 

Mi = v ’ «y(r Q , /3) 

7TZ 2 


( 4 ) 


where 
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^(r 0 ,/3) = 


cos| — — 

,2 R 


\ 


exp 


-X 1 


'/r^ + r^ - 2rr cos /3 + z^ / 
o o / 


[n 


r dr d/3 
o o 


( 5 ) 


The neutral atom density, n Q , at the point (r,z) is calculated from 
the neutral arrival rate by, 

M 0 ( r > z ) 

n (r,z) = (6) 

v o 

This relation assumes all trajectories are parallel. It will yield 
slightly lower neutral density estimates close to the thruster where the 
applicable relation for /j, Q = i n Q v o ; v Q is the thermal velocity, 
y(8kT/j7m . 

The neutral atom arrival rate is determined by the relation 


|U 0 (r,z) = - f(x)] (7) 


where 

C ^ fraction of average ion emission rate that represents neutral 
emission rate 

— 2 - ~ ~ 

v~ ~ " average ion emission rate" T~/'q7fRg 

Fjg view factor from thruster to point (r, z) 

f(x) fraction of total neutrals emitted which have undergone charge 
exchange up to axial station, x 
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The view factor Fjg for a disk source, assuming uniform neutral 
emission rate across the surface is, (ref. 13) 



where x = Z/R q , and y = r/z. 

Substituting (4), (6), and (7) into (2) and indicating integration over 
all variables (r Q from o to R q ; 0 from o to 2jr; r from o to <»; 
z from o to z) yield 


r 

* 

J 

-2. 16C(n,x) 

. « 2 - . 


C 


(1 - too ). /he 

v 8kT 


2nr dr dz 


(9) 


But, the fraction of total neutrals emitted which have reacted up to dis- 
tance x is 


N(x) 

f(x) = ? 

< C l" 


( 10 ) 


Combining (10) and (9) and converting to dimensionless variables yield 


r 

r 

f 12 (i -f(x)) 

'X=0 

Jy= 0 

A 


y dy dx 


(ID 


where 
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A = 


8. 64 a 10 C(n,A)i/R o 



I(w,/3) 


^)(r 0 ,/3) 



( 12 ) 


w 



Equation (11) was solved numerically for various values of A. 
Since A (eq. (12)) contains several variables it was evaluated for one 
set of conditions. It can then be expressed in terms of ratios to these 
” standard” values. 

The following values might be typical of a 30 -centimeter diameter 
thruster using mercury propellant: 

R q = 15 cm T = 500 K C(n,A) = 10.92 

I =1.5 a m = 200 (Hg) a 10 = 6x10” 15 cm 2 

The resulting value of A for these conditions is 1. 56. Thus, 


n = 1 
A = 25 


- a = i. J - ai ° m 

\6X10“ 15 / ^ 10 - 92 / V 0 - 00212 / ' T ' ^ 20o/ V 15 / 

Charge -exchange cross-section values are available for some of the 
reactions of interest (refs. 14 to 17), but have only been estimated for 
others of primary concern (ref. 18). 
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A plot of the fraction of emitted neutrals which have engaged in 
charge exchange as a function of the distance from the thruster is shown 
in figure 1. The corresponding charge exchange formation rate at any 
location is shown in figure 2. In these plots several A -values are 
shown. A change in A-value can represent a change in any one or sev- 
eral of the parameters indicated in equation (13). Several X -values are 
also shown, representing the range from a fairly collimated beam 
(n = 1, X = 60) to one more divergent (n = 1, X = 25). A caution in inter- 
pretation is to be noted here. While it appears that a less well colli- 
mated beam has more charge exchange than a more directed beam (com- 
paring, say curves for A = 1 for the (n = 1, X = 25) and (n = 1, X = 60) 
conditions) it should be remembered that A contains the parameter 
C(n,X). The same value of A for these two cases, then, does not 
represent otherwise identical conditions. When this is taken into ac- 
count, it will be found that there is little difference in charge exchange 
formation in comparing just the effect of beam divergence. Values of 
C(n,X) may be found in reference 9. 

Some other points that may be noted in figures 1 and 2 are: 

(1) Most of the charge -exchange interaction that occurs will have 
taken place by about 8 to 10 thruster radii downstream. 

(2) Beyond about 8 thruster radii, the charge exchange formation 
rate varies inversely with the square of the distance (fig. 2). 

(3) Below about 10 percent reacted, the charge -exchange formation 
is proportional to A. Fraction reacted at some value of A not plotted 
is thus readily estimated 


f(x) - [fMj^ 
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CONCLUDING REMARKS 

Each spacecraft layout involving electrostatic thrusters presents 
its own unique set of conditions which should be systematically checked 
for possible interaction between the ion beam and spacecraft compo- 
nents. For this reason, this report has presented an outline of the 
procedure that has been followed by the author at the Lewis Research 
Center to check the possible spacecraft contamination from ion thruster 
operation. 

Details of the calculation procedure have not been presented if they 
were available in already published references. The references have 
been cited in each instance . 

Details and graphical results of charge -exchange calculations were 
presented since this particular material has not been published 
previously. 
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APPENDIX - SYMBOLS 

constant, defined by eq. (12) 
constant in ion profile equation 

fraction of average ion emission rate that represents neutral 
emission rate 

view (or configuration) factor 

fraction of emitted neutrals that have exchanged charge 

total current, a 

integral in eq. (11) 

integral in eq. (4), defined in eq. (5) 

current density on center line, e = o 

current density in radial plane at angle e from beam 
centerline 

Boltzmann constant 

particle mass 

total number of ions exchanging charge 
constant in ion flux eq. (1) 
neutral particle density 
electron charge 
radius of thruster 

radial distance in plane downstream of thruster 

radial distance in thruster exit plane 

area 

temperature, K 
volume 



particle velocity 

dimensionless radial variable = (r /R ) 
dimensionless length variable = Z/R Q 
dimensionless radial variable = r/Z 
axial distance from thruster exit 
angle variable (sketch (a)) 

angle variable in thruster exit plane (sketch (a)) 

angle between thruster centerline and trajectory leaving thruster 
exit center 

constant in ion beam profile eq. (1) 
ion arrival rate per unit area 

average ion emission rate from thruster = I a /q7rR Q 
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Figure 1. - Variation of charge exchange neutral 
fraction with distance from thruster exit. 

(n, X) determined by equation (l); A by equation 
( 13 ). 
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Figure 2. - Variation of charge-exchange formation 
rate with distance from thruster exit, (n, l) 
determined by equation (l)j A by equation (13). 
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